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Abstract. 



> 

ff^ ' We consider a simple physical model for an evolving horizon that is strongly 

op . interacting with its environment, exchanging arbitrarily large quantities of matter with 

its environment in the form of both infalling material and outgoing Hawking radiation. 
We permit fluxes of both lightlike and timelike particles to cross the horizon, and ask 
, how the horizon grows and shrinks in response to such flows. We place a premium on 

providing a clear and straightforward exposition with simple formulae. 



o 
o 

^ 

o 

O ■ To be able to handle such a highly dynamical situation in a simple manner we make 

I , one signiflcant physical restriction — that of spherical symmetry — and two technical 

^Jrj' mathematical restrictions: (1) We choose to slice the spacetime in such a way that the 

space-time foliations (and hence the horizons) are always spherically symmetric. (2) 
Furthermore we adopt Painleve-Gullstrand coordinates (which are well suited to the 
' problem because they are nonsingular at the horizon) in order to simplify the relevant 

calculations. Of course physics results are ultimately independent of the choice of 
coordinates, but this particular coordinate system yields a clean physical interpretation 
of the relevant physics. 

We find particularly simple forms for surface gravity, and for the first and second law 
of black hole thermodynamics, in this general evolving horizon situation. Furthermore 
we relate our results to Hawking's apparent horizon, Ashtekar et aZ.'s isolated and 
dynamical horizons, and Hayward's trapping horizon. The evolving black hole model 
discussed here will be of interest, both from an astrophysical viewpoint in terms of 
discussing growing black holes, and from a purely theoretical viewpoint in discussing 
black hole evaporation via Hawking radiation. 
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1. Introduction 

Throughout the 1970's the black hole spacetimes considered in the scientific literature 
were typically globally stationary and asymptotically flat — with true mathematically 
precise event horizons [absolute horizons]. When accretion was considered it was most 
often handled in the test-matter limit, where the self-gravitation of the accreting matter 
was ignored. After the accretion rate was calculated in this approximation one might 
feed this back into the black hole geometry to allow the spacetime to evolve in a 
quasi-static manner, but the horizon itself was almost always treated as though the 
idealized Schwarzschild or Kerr event horizon captured almost all of the relevant physics. 
The "Membrane paradigm" is perhaps the apex of this quasi-stationary approach [1]. 
Similarly, calculations concerning Hawking radiation are typically carried out in the 
test-field limit, where the gravitational field generated by the outgoing Hawking fiux is 
ignored [2] . Once the Hawking fiux is calculated, this is set equal to the mass loss rate 
of the underlying black hole, with the black hole again evolving in a quasi-stationary 
manner [3]. 

Those explicit models where exact solutions of the Einstein equations were used 
to investigate black hole horizons beyond the quasi-stationary approximation generally 
made rather strong assumptions on the nature of the matter crossing the horizon. For 
instance, in Oppenhcimer-Snyder collapse [4] one is limited to dealing with zero pressure 
dust, while in the Vaidya solutions [5] one is limited to lightlike radiation crossing the 
horizon [6] (either incoming radiation, or outgoing radiation, but not both). 

More recently, considerable work has been done on relaxing these assumptions in 
various ways. Motivated partly by advances in numerical relativity, Ashtekar et al. [7], 
and independently Hayward [8], have developed formalisms that address what they 
variously call "isolated", "dynamical", and "trapping" horizons. These modifications 
of the notion of event horizon exhibit some but not all of the properties of the 
"apparent" horizon defined (for instance) in [9]. (For more technical details, see the 
recent review [10].) Herein we shall present a particularly simple framework that, while 
assuming spherical symmetry, does not require asymptotic fiatness (or more importantly 
any notion of stationarity) either of the global spacetime or of the black hole itself. 
Similarly, no a priori constraints are placed on the nature or quantity of matter crossing 
the horizon, in either direction. Thus the framework provides a useful illustration 
and testing ground for some of the ideas recently popularized by Ashtekar et al. and 
Hayward. 

The choice of Painleve-Gullstrand coordinates (t, r, 9, (f)) rather than the more 
familiar Schwarzschild coordinates (t, r, 9^ (p) is motivated by a desire to clarify issues 
such as how the black hole is able to grow despite the supposed "frozen" nature of 
Schwarzschild time at the horizon, and the related issue of wanting to calculate quantities 
at the horizon using a coordinate system that does not break down at the horizon. The 
use of t distinguishes the Painleve-Gullstrand time coordinate from the Schwarzschild 
time coordinate t. These two coordinates agree at spatial infinity where they correspond 
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to the proper time of observers at infinity. We do not, however, require that such an 
asymptotically flat region need exist. 

We also include some discussion of the relation between our evolving horizons, 
apparent horizons, and the isolated/dynamical/trapping horizons of Ashtekar et al. and 
Hayward. It has recently been conjectured that in a fully quantum gravitational universe 
true mathematically precise event horizons will not form, only long-lived apparent 
horizons [11]. This fact has led to claims of a solution to the black hole information 
paradox [11, 12, 13, 14]. While the present work does not directly shed any light on this 
issue, it does provide a classical arena wherein such ideas can be discussed. | 

The organization of the paper is as follows. In section 2 we derive the form of the 
time dependent, spherically symmetric metric in Painleve-GuUstrand coordinates and 
in sections 3 and 4 we use this formulation to derive equations for the rate of change of 
the horizon radius and area. In sections 6 and 7 wc calculate various fluxes and stresses 
at the horizon using a particular choice of tetrad basis. In sections 5 and 8 we relate this 
work to that done by Ashtekar et al., and Hayward, and discuss the unified formalism 
for expanding, stationary, and evaporating horizons. 



2. Painleve-GuUstrand coordinates 



It is well known that any spherically symmetric gravitational field can always be written: 

ds^ = g^^{r^ t) dt^ + 2grt{r, t) drdt + grr{r, t) dr^ + R{r, tf d^f . (1) 

The remaining coordinate freedom in the r-t plane can then be used to reduce the four 
independent functions above to two. For instance, it is also well known that a static, 
spherically symmetric gravitational field can always be written, at least locally, in terms 
of curvature coordinates [aka Schwarzschild coordinates), as [16, 17] 



ds' 



_g-2$(r-) 



1 - 



2m (r) 



dr^ 



+ dQ^ 



(2) 



1 — 2m(r)/r 

and furthermore for static spacctimes this representation is quite useful — and in the 
absence of horizons is even globally useful. Then for a time- dependent spherically 
symmetric gravitational field it is natural to write 



ds' 



_g-2$(r-,t) 



1 - 



2m(r, t) 



dt' + 



dr^ 



1 — 2m(r, t)/r 



+ r' dQ^ 



(3) 



Indeed, we know that any general spherically symmetric metric can be written in terms 
of two unknown functions and here we are just making them time dependent. So no one 
can stop you from adopting such a coordinate system, it's just that (as we shall soon 
argue) this particular coordinate system is less useful than one might at first suppose. 

Now for any given slicing of spherically symmetric spacetime into hypersurfaces, and 
in particular for the choice above, one can look for the existence of apparent horizons 

I For a related quantum framework phrased in terms of the spherically symmetric sub-sector of the 
coupled Einstein-scalar field system see [15]. 
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in the standard way (see, for instance, [9]). That the apparent horizons depend on 
the shcing of spacetime was shown most clearly in reference [18] where non-symmetric 
slicings of Schwarzschild spacetime were found for which no apparent horizons exist. 
An explicit example of this can be found in [19]. Because of this particular "feature" 
of apparent horizons both Ashtekar et ai, and Hay ward, when discussing general 
spacetimes, have to introduce considerable technical machinery to avoid pathologies 
due to potentially ill-behaved slicing of the spacetime. 

In contrast, in this article we will restrict ourselves to spherically symmetric 
spacetimes and spherically symmetric slicings. So much, though not all, of the technical 
superstructure simplifies. (The whole point of the present article is that we want 
to retain just enough technical machinery to do the job, and eliminate as many 
complications as possible.) Then for the coordinate choice in equation (3) the apparent 
horizon is informally defined by the implicit relation 2m{r,t)/r = 1. (We will be more 
precise and formal below.) That is: 

2m{rH{t),t) =rH{t). (4) 

In the Schwarzschild metric m{r,t) = M, where M is a constant, and this just 
corresponds to the familiar Schwarzschild radius. But once you start asking questions 
about the evolution r//(t) of this apparent horizon you are plagued with multiple divide 
by zero errors, as the matrix of metric coefficients, when written in these particular 
coordinates, is a singular matrix at the horizon. The ultimate reason for this behaviour 
is that Schwarzschild coordinates are reasonably good (or at least not intolerably bad) 
for probing the geometry of a static black hole, but are not particularly good (in fact, 
downright awful) for probing the future horizon of an evolving black hole. 

To deal with a dynamical black hole, we should use a coordinate system that 
is well behaved at the apparent horizon. For this article we will use Painleve- 
Gullstrand coordinates which are particularly simple and lead to a particularly nice 
physical picture. (Other non-singular coordinate choices such as Eddington-Finkelstein, 
generalised Vaidya, or double-null coordinates are possible and the physics is entirely 
equivalent. § We concentrate on Painleve-GuUstrand coordinates because they give a 
particularly nice form to the equations.) 

Starting from (3) we perform the change of coordinates t — > i{t, r) where i will be 
the Painleve-GuUstrand time. Thus 
di dt 

dt = — dt + — dr = tdt + i' dr. (5) 
ot or 

Substituting in for dt. 



d.^ = -e--(-*) ( 1 - '-^] (idi-tdr]\ + r^dn^ 

\ r J \i t J l-2m(r,t)/r 

Expanding out and demanding that = 1 gives the condition 



(6) 



1 - 2m{r,t)/r ' ^ ' 

For an early article presenting an interesting formulation in terms of double null coordinates see [20] . 
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That this partial differential equation, which is simply a first-order linear homogeneous 
equation for t{t, x), always has a unique solution can be shown [for instance] by the 
method of characteristic curves [21] ||. Imposing (7) now leads to 

p-Hr,t) 



9tr 



t 



-\/2m{r^ t)/r, 



which can be written as 

9tr = ±c(r, i)^2m{r,i)/r = v{r, i). 
This imphcitly defines c(r, i) and leads to 



(8) 



(9) 



gu = -ic'-v'). (10) 

Thus the final metric, now in Painleve-Gullstrand form, becomes 

ds2 = -c(r, ty dp + [dr + v{r, i) di]^ + r^ dQ\ (11) 
or equivalently 

ds^ = -[c(r, i)^ - v(r, tf]dP + 2v(r, i) dr dt + dr^ + r^ dQ^. (12) 

That is, any spherically symmetric space-time, regardless of whether it is static or not, 
can always locally be put into this form. Note that, as is usual in Painleve-Gullstrand 
coordinates, surfaces of constant t are spatially fiat, and there is a natural notion of 
"outwards" and "inwards" associated with increasing or decreasing the r co-ordinate. 
Furthermore, four- dimensional asymptotic flatness would correspond to imposing c — 1 
and f ^ as r ^ oo. (We will not need to impose asymptotic flatness in the discussion 
that follows.) Note also that the static Schwarzschild case would be given by setting 
c — 1 everywhere, while v{r) = y^M/V, with M a constant. Explicitly we have 



and 



9ab — 











Vi 




hij 



















(13) 



(14) 



where hij is the metric of flat Euclidean 3-space in spherical polar coordinates. The 
ingoing and outgoing radial null curves are defined by ds^ — and equation (12) gives 
dr 

= -v{r,t)±c{r,t). (15) 

dt 

Thus we can define the location of what we will call the evolving horizon by the very 
simple and intuitive condition 

c(r,t) = ?;(r,t), (16) 

II Of course this uniqueness holds only subject to the specific choice grr = 1, and a suitable boundary 
condition (such as t' ^ t when r — > oo). What is certainly not unique is the general process of 
constructing a nonsingular coordinate system at the horizon. Other coordinate systems, such as 
Eddington-Finkclstcin coordinates, are possible and will lead to similar results, that are likely however 
to differ in small technical details. 
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which has the simple physical interpretation that when c{r,i) < v{r,i) the outgoing 
light ray is being dragged backwards to smaller values of r. This implicitly defines a 
function rH{i) such that 

v{rHit),t)^c{rH{i),i), (17) 

which is equivalent to 

2m{rH{i),i)=rH{i). (18) 

Now the question is; what is the evolution of this function rnii) in terms of the stress- 
energy at rH{i)l 

3. Surface gravity and the First law 

Using Maple (or some equivalent) to compute the Riemann tensor of the metric (12) in 
orthonormal components we calculate 

R (19) 



so that m{r,t) can be physically identified as the Hernandez-Misner mass function [22]. 
Similarly, since we are dealing with a spherically symmetric spacetime, we could compute 
the Hawking-Israel quasi-local mass function 

mui{r,t) = ^ [l-g'^'VarVbr] , (20) 

which again leads to the mathematical quantity m(r, t) appearing in the metric being 
physically identified as the "mass inside radius r at time t" . 
Now since 2m(r, t) = r at the evolving horizon, we have 

2m{rH{i),i) + 2m'{rH{i),i)fH{i) = fH{i), (21) 
which we can recast as 

m{rH{t),t) = ^ )_E\1lJ1. r^it). (22) 

Introducing Ah = 47rr|^, the area of the evolving horizon, we have 
. / /~x -A 1 [l-2m'(rH(i),i)] : ,~, 

where we have seen that the function m(r) can be physically interpreted, due to the 
spherical symmetry, as the mass contained within a radius r. Therefore this has exactly 
the form of the first law of black hole mechanics, dm = -^k, dA, though now in a 
completely dynamical context, provided we agree to focus on partial derivatives with 
respect to i ^, and also to identify 

^H{t) = (24) 

^ If we were to compute the total time derivative of m{rH{i),i), we would obtain 

dr 2, 
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as the "surface gravity" . + To justify this interpretation of kh, it is useful to first define 
the outward radial null vector 



(l,c{r,i)-v{r,i),0,0) 



c{r,i) 



(25) 



and verify that gab ^" — 0. The overall normalization is chosen to make things simpler 
below. It is also useful to define the inward radial null vector 



n — 



l,-c{r,i)-v{r,i),0,0 



0, and that gab = —2. (This choice of 



c(r, t) 

and verify that this is also null gab n"' 
normalization is the most "symmetric" we have been able to find.) 
Because of spherical symmetry we must have 

t Wbt = Ke T; V^n" = n"; 



(26) 



(27) 



where the scalars and k„ are defined everywhere on the spacetime, not just on the 
evolving horizon. Computing and k,„ for the given and yields 

c'(r, t) — v'{r, t) 



Ke{r, t) 



c(r, i) 

1 [2m{r,i)/r -2m'{r,i)] 



2r 



+ 



c'{r,t) 
c{r, i) 



Kn{r, i) = 



•\/2m(r, i)/r 

—c'{r, i) — v'{r, i) 
c(r, i) 

1 [2m{r,i)/r -2m'{r,i)] ^ c'(r,t) 

c(r, i) 



1 - 



2m (r, t) 



2r V'2^^^M)7^ 

On the evolving horizon Ke reduces to 

~ l-2m'{rH{i),i) 
Knit) = - 



+ 



1 + 



2m(r, t) 



(28) 

(29) 

(30) 
(31) 

(32) 



2r^(t) 

a formula which has all the appropriate limits and is compatible up to normalization with 
the results in reference [17]. In particular, if the geometry is static (time-independent) 
then it is obvious that evolving, apparent, event, isolated, dynamic, and trapping 
horizons coincide, and that kh as defined above reduces to the standard definition. 
(In reference [17] the spacetime was in addition assumed asymptotically fiat, which we 
do not need to assume in the current analysis. The asymptotic fiatness was then used to 

This result, though simple, does not appear to be particularly useful, and does not seem to have a 
straightforward interpretation in terms of black hole thermodynamics. The best one can apparently do 
in this regard is to rearrange the above as 

dm{rH{i),i) _ rh{rH{i), i) 

di ~ 1- 2m'{rH{t),i)' 

+ For a discussion of situations wherein the "entropy= area/4" law might fail, see for instance [23, 24]. 
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motivate a particular non-local normalization for surface gravity. The current definition 
of "surface gravity" is in contrast local, and only uses information that can be extracted 
from the vicinity of the evolving horizon itself.) 

Now we want to be dealing with the outermost evolving horizon, which for simplicity 
we define in the purely intuitive sense. Thus for purely kinematic reasons 

^^^^<1 for r>rH{i), (33) 
r 

which implies 

/2m 



<0 ^ l-2m'{rH{i),i)>0. (34) 

H 

Note that in the case of a spacelike evolving horizon the horizon might, in principle, 
intersect the constant time slices many times. In such a situation this condition will 
only identify the outermost intersection. See the appendix on the Vaidya spacetime for 
more details on this point. Apart from the exceptional case 1 — 2m'(rj^(t), t) = 0, we see 
that the surface gravity on the outermost horizon is guaranteed to be positive. Indeed, 
from the above we can see that the case where 1 — 2m'{rH{t),t) = corresponds to 
i^H{t) = 0, the condition for an extremal horizon.* For a discussion of the subtleties 
that can occur once one has multiple nested dynamical horizons, see [25]. 

Aside: We could also calculate k„, corresponding to the ingoing null geodesies, on the 
evolving horizon. Though there is no difficulty in doing so, there is no clear physical 
interpretation of the resulting quantity: 

Mr^(t),t) = «:^(t) + ^^^^:^. (35) 

c{rH{t),t) 

Aside: From the defining relations for K£ and Kn, the spherical symmetry of the 
spacetime, and the normalization relation between i and n, it is easy to see that 

V^r = -Knt] t VbU^ = -Ken''] (36) 

or equivalently 

i^n = ^ ; i^e = 2 ■ ^^^^ 

Aside: From the definition of K£ 

t Vbt - t (38) 
we see that rescaling the null vector i ^ ai (which implies n a^^n), will result in 

^^{at) = a + i ■ Va. (39) 

* For example, in the case of the Reissner-Nordstrom spacetime we have m{r) = M — Q^/2r and a 
Uttle bit of algebra is enough to show that the condition 1 — 2m' [rn) = is equivalent to the condition 
= Q^, the usual condition for an extremal Reissner-Nordstrom black hole. However in more general 
situations the vanishing of the surface gravity is typically taken to be the primary definition of what 
one means by extremality. 
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So at the evolving horizon, where £ • V — > dt, we have 

CK 

KH{a i) =a Knit) + — • (40) 

Ch 

That is, among all possible normalizations for (and so implicitly for n), the one we 
have chosen is seen to not only be "symmetric" but additionally to eliminate any explicit 
time derivatives at the evolving horizon. 



4. Null energy condition 

Now calculate the quantity Gab ■ Using Maple (or some equivalent) it is easy to see 
that 



c(r, t) H y/2m{r, t) /r c{r, t) r 
which we can rearrange to yield 



m(r, i) = ^c(r, i) T t - c'{r, i) ^2m(r, i) r (^1 - ^J^^j ■ (42) 

This mass formula applies at any value of r. In particular at the evolving horizon we 
have the very simple result 

rh{rH{i),i) = ^rl{i) c{rH{i),i) Gab 1 1. (43) 

Invoking the Einstein equations Gab — StJ" T^ab (with Gjv ^ 1), we see 

m(r^(f), f) = 47rr|(f) c(r^j(f), t) Tab t (44) 

Note that this only includes the mass change due to flux across the instantaneous 
location of the evolving horizon, this does not yet include the effect due to the motion 
of the evolving horizon. 

But in view of equation (22) we have 

. - ^ 2m{rH{i),i) ^ STrrlji) c{rH{i),i) Tab ^ i' 

1 - 2m'(rH(t),t) 1 - 2m'(rH(t),t) ' ^ ' 

Furthermore since we have already seen that 1 — 2m' {r nit) it) > for any non-extremal 
outermost horizon, we deduce that the sign of fn is the same as the sign of Tab But 
this is exactly the combination that enters into the Null Energy Condition {Tab -^^ ^ 0; 
the NEC) evaluated on the apparent horizon. 

That is, as long as the NEC is satisfied the horizon cannot shrink. This is compatible 
of course with the standard analysis in terms of the Raychaudhuri equation [9], but 
here we can sec the result dropping directly out of one component of the Einstein 
equations applied to the evolving horizon, (i We emphasise that the NEC, because it 

)j Similar violations of the NEC occur at the throats of dynamically evolving wormholes [26] — and 
one also encounters problems similar to those for evolving horizons if one tries to describe evolving 
wormholes in Schwarzschild curvature coordinates. Schwarzschild curvature coordinates are ill-adapted 
to discussing the evolution of a wormhole throat. 
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is the weakest of the standard energy conditions, leads to the strongest form of the 
singularity theorem [27, 28]. Because the singularity and area increase theorems can 
both be phrased in terms of the NEC we do not really need to consider the WEC, SEC, 
or DEC. 

Aside: We can also define the total quasi-local mass inside the evolving horizon by 

mnii) = m{rH{i),i) = r^f(t)/2, (46) 

and then 

dmff ^ m{rH{i),t) ^ ^Trrjjji) c{rH{i),i) Tg^ 1 1 
dt 1 - 2m'(rfl(t),f7) \ - 2m' {r n(i) ~t) ' 

This evolution law now characterizes the total mass change, now with contributions 
both from mass fiux across the evolving horizon and from the motion of the evolving 
horizon itself. 

Aside: In contrast, consider the ingoing radial null direction. For G^h n"" n'' we have 
the result 

c{r,t) r^^2m{r,t)/r c{r,t) r i v „ ; 

which we can rearrange to yield 



m(r, i) = ic(r, t) r' ^ n« n' - c'{r, i) ^2m(r, i)r(^l + ^ j , (49) 
but this is nowhere nearly as useful since its form on the evolving horizon is not as nice. 



5. Apparent, Isolated, Dynamical, Trapping, and Evolving Horizons 



While we have used Painleve-GuUstrand coordinates above to get some very simple 
and explicit results, we note that similar but coordinate independent definitions of 
various types of horizon are given by [7, 8, 9]. To sec the relation between the 
evolving horizons at 2m(r) = r and the apparent /isolatcd/dynamical/trapping horizons 
of Hawking/Ashtekar/Hayward we need to calculate the expansions 6i and 6n of the 
radial null vectors and n". Defining the expansion of the outgoing and ingoing radial 
null curves as 



9"" + 



,ab 



v„4 = v„r - K,- 



VaTlfe = Van" - Kn, 



(50) 



(51) 



we get 
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rc r 




(53) 



So we see that 9£ changes sign exactly at the evolving horizon rH- This is enough to 
see that the evolving horizons of this article coincide with the standard definition of 
apparent horizon — indeed all you really need for this to hold is to have a spherically 
symmetric spacetime with a time slicing that respects the spherical symmetry. 



Now define an "Ashtekar horizon" H as follows: 

i. H is a three-dimensional, timelike/ null/ spacelike hypersurface. 

ii. On H the expansion of is negative ^„ < 0. 
Hi. On H the expansion of is zero 6^ — 0. 

The null case represents an "isolated horizon" and the spacehke case represents a 
dynamical horizon. (This is the astrophysically relevant case where accretion dominates 
over Hawking evaporation.) The timelike case corresponds to a situation where Hawking 
evaporation dominates over accretion — this might be referred to as an "evaporating 
horizon". From the explicit formulae for the expansions 6£ and 6n it is clear that our 
evolving horizon is also a horizon in the sense of Ashtekar et al. 

Finally, a "Hayward horizon", or more precisely a future outer trapping horizon, 
replaces condition i with 



CJf < 0. 



But explicitly 



CJi = (n • V) 



and so on the evolving horizon 



1 - 




(n- V) 



2m 
r 



H 



- — (n- V) 
rH 



2m 
r 



H 



On the evolving horizon we also have 



whence 



That is 



(n- V) 



2dr 



1 



2m f2m\' 

2 — 

cr \ f J 



H 



' H 



+ 2 1 

c 



2m') 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



H 



But we have already seen, on purely kinematic grounds, how to relate m\H to fn, so 
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We have furthermore already argued that 1 — 2m' \h > 0, so under the rather mild 
condition that rn > —2ch we have CnOi\jj < 0, in which case our evolving horizon is 
also a Hayward-style future outer trapping horizon. 

In the situation where rfj < —2ch we have a somewhat unusual type of spacelike 
evolving horizon — in the present context this would correspond to rapid "evaporation" 
of the horizon. Since the evolution of the horizon now lies outside the hghtcone, 
evolving forward along curves of n" does not take us from a region of positive $£ to 
negative O^. Timelike observers are still instantaneously forced to move inwards, but 
they immediately escape the horizon — so that subsequently some of them may turn 
around and escape to "infinity" . Indeed, once the observer exits they will not be able 
to re-enter the horizon unless the rate of evaporation slows down. 

In short, the very simple and intuitively clear definition of evolving horizon that 
wc advocate in this article is compatible with all the standard definitions (apparent, 
isolated, dynamic, evaporating, trapping) common in the literature. The one thing the 
evolving horizon is not, is that it is definitely not an event horizon [absolute horizon]. 
And one very important message to take from this entire discussion is that event horizons 
[in the strict technical sense] are not essential to investigating black hole physics [11, 29]. 



6. Energy momentum in a natural tetrad basis 

To understand the energy-momentum fiuxes and stresses at the evolving horizon it is 
extremely useful to adopt a suitable orthonormal tetrad basis. Once you have picked a 
tetrad, which is effectively a choice of privileged observer, you can begin to ask questions 
about densities, fluxes, and stresses measured by that observer. 

To construct an appropriate "natural" tetrad, start with the observation that the 
radial null curves are given by 

%^-v±c. (61) 
which implies that the vector (1, —v. 0, 0)°^ is certainly timelike. In fact 

^„ ^ (1,-^,0,0) ^ + ^g2) 

is easily seen to be a timelike unit vector. (In an analogue interpretation of the metric, V 
would be the 4- velocity of the "medium" , and our natural observers would be co- moving 
with the medium [30].) Similarly, 

5« = (0,1,0,0); Sa = iv, 1,0,0); (63) 

is easily seen to be a spacelike unit vector that is orthogonal to V. We take these as the 
first two elements of the tetrad. Specifically we set 

ea"= (n-^",^",*^")- (64) 

Note that 

V- + S^^ ~ °' °^ = r ; (65) 
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(66) 



where i is the outgoing null vector we previously chose to make the surface gravity 
computation simple. One may easily compute 



(y-V)V^-S; (V-V)S^-V; 

c c 



(67) 



(68) 



Thus an observer who "moves with the medium" is not in geodesic motion except in 
the special case d — 0. 

Writing = and r" = S"- and using this natural tetrad (and its inverse as 
determined by Maple jf), the orthonormal components of the Einstein tensor (at any 
point in the spacetime) are easily calculated to be 



_ 2m' 2c' 2m 

^it — — T ; 



G 



tf 



2d 2m _ 
cr \ r ^ 

2m' 2d 2m 
+ — + 



' cr ' cr^ ^2m/r 



The only "complicated" component of the Einstein tensor is the transverse one: 

+ - 1 + - 

cr cr \ r / 



Gaa — Gix — 



m c 
r c 



2m ^ 



(69) 



(70) 



(71) 



(72) 



1 2m ( m 
" \ Amcr 



1 + 



2dr 



rm 



m 



1 



m' 2d 2cd 

— + - 

m c c^ 



The Einstein tensor at the evolving horizon simplifies rather drastically 



2m' 2d 

^ + — ; 



G 



if\H 



G, 



ff\H 



cr 
2d_ 

1 

cr 

2m' 2d 2m 
+ — + 



cr cr^ ' 



(73) 



(74) 



(75) 



ft Wc particularly wish to warn readers against over-enthusiastic use of the frame command in Maple, 
or its equivalent in other packages. The frame function will provide some orthonormal tetrad for the 
specified metric, but the orthonormal tetrad it provides is certainly not unique, and may not always 
be the most useful. In particular, if the spacetime contains any form of horizon, the vectors in the 
orthonormal tetrad provided by the frame command qiiite often have divergent components at the 
horizon. The fact that the tetrad is not unique is of course an unavoidable consequence of local Lorentz 
invariance, and the divergence of the components of the tetrad is a side effect of attempting a "t; = c" 
Lorentz transformation at the horizon, so similar problems will show up in any symbolic computation 
package. We have found it best to determine a suitable well-behaved tetrad "by hand", and then 
explicitly feed it to the symbolic manipulation program for further computations. 
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with the only "comphcated" component being 

m" 3c'm' _|_ 3c' ^ m 
^ r cr 2cr 2cr^ 



G 



IH 2m' 



1 

+ - 

c 



m ^ c cc 



(76) 



Now write the Einstein equations as Gab 
Then (at the evolving horizon) 



PH 



in 



Pr,H 



-tr 



J. T'i 



Git ^ 
Stt 

Gjf _ 
Stt ~ 



m 

c' 



+ 



47rrc' 



47rrc 
m' 



+ 



Stt 47rr^ ' 47rrc 
with again the only real complication being 

G,-^ m" 3c'm' 3c' 
- Jflfl - J-si - -TT- - - H 



c 

SttGtv 2a65 and adopt units where Gat = 1. 



(77) 
(78) 
(79) 



+ 



m 



Stt 



+ 



Stt r 

m 



Stt cr IGvr cr 

. 2c'r ^ ; 
1 + 2m' 



(80) 



+ 



87r c 



m 
r 



^ c cc 



Note that these must be interpreted as the density, flux, and pressures as measured by 
the particular observer moving with 4-velocity V"-. In terms of these quantities we have 
the simple results 

m = 47rr'c {pn + Pr,H + ^fn) = ^Tir^c [Tabti%, (81) 

and 

m' = Anr^c [pn + fu) = ^i^r^c [T,,rV%. (82) 
This now implies that the evolving horizon shifts according to the formula 
^T:cHrl{pH + Pr,H + V h) Sncnrjj [Tabt^U 



and similarly 



rriH 



1 - %T^rl{pH + fH) 1 - 8nrjj[TabV-i%' 

AnCHrjj{pH+Pr,H + 2fH) ^ ATTCHrjj [Tabi''i% 

l-87rrl[T,bVH%- 



(83) 



(84) 



l-S7CrUPH + fH) 

Furthermore the surface gravity is 

i-w i-sMph + Ih) 1 - 8xir„ii/°<»]„ 

This now provides explicit formulae for the evolution and properties of the horizon in 
terms of various quantities (densities, fluxes, pressures, etc..) evaluated at the horizon. 
The only real deflciency in the tetrad we have chosen (and this is a matter of taste, not 
a matter of physics), is that while it is easy to write down, and easy to interpret, and 
easy to calculate with, the 4-velocity V is not geodesic (except in the rather special case 
where c{r,t) ^ c is a constant). Now in that special case not only is V geodesic, but 
also / = 0, the net flux seen by a freely falling observer located at the evolving horizon is 
zero. This leads us to ask whether this behaviour can be emulated by a slightly different 
choice of tetrad? 
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Consider now the 4- velocity 

ya 1 



[l,-(^; + c/3),0,0], 



(86) 



where |/3(r, t)| < 1. This corresponds to an observer who is infalhng with 3-velocity 

vs{r, t) = v{r, t) + c(r, t) /3{r, t). (87) 

It is easy to check that y is a unit timehke vector, and easy to construct an appropriate 
spacehke vector that is orthonormal to it: 

1 



5« = 



[-f3,c + vP, 0,0] 



(88) 



This now allows us to construct a generic non-null tetrad (one element of the tetrad, V, 
is guaranteed to be timehke, the other three spacehke) 



In terms of this tetrad one now has 



p = 



87r 



1 



Pr 



Anrc{l - /32) 
87r 

1 

47rrc(l - /?2) 

Gca Qb 



2m 

r 



2m 
r 



-13 dt 



2m 

r 



/5 + 



2m 



Sn 



1 



m 



2m 



c{l-(3')—+\l+(3J— \ c' + dt\J — 



47rrc(l - /?2) 

Then at the evolving horizon, where 2m = r, we have 



2m 



Ph 



fH 



P 



H 



1 


47rrc(l — 


/32) 


1 




47rrc(l — 




1 




47rrc(l — 





m 



c{l-(3'')— + {l + (3fc' + 



,m 



H 



m 



-/3--(l + /3)'c' 
r 

-c(l-/3^)^ + (l + /3)'c' + 



H 



m 
r 



(89) 

(90) 
(91) 

(92) 
(93) 

(94) 
(95) 

(96) 
(97) 

(98) 



H 
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The special case /3 — > reproduces the results based on our "natural" tetrad. From 
these results we can see that the notion of "flux across the horizon" is both subtle and 
less useful than one might at first suppose — the flux depends explicitly on the tetrad 
(effectively, one's choice of privileged observer), and indeed the "flux across the horizon" 
can in many situations be made to vanish by choosing a suitable observer: 



/3h 



[1 + Ph] 



rc 



m 



(99) 

H 

But even though the notion of flux is somewhat more subtle than expected, the evolution 
of the horizon as encoded in quantities such as m|ij and drn/dt is perfectly well defined 
and easy to deal with. 

8. Formalism for generic horizons 

We now show how to write the change in the area of the horizon as a function of the 
expansions defined above. First, choose a coordinate basis {6, 0) for the surface of the 
horizon at a given instant in time (on a given slicing). Then the area of this 2-surface 
will be 



Ah= \/det h de d(p, (100) 



where h is the induced metric on the 2-sphere given by its embedding in the full 4D 
spacetime. For a spherically symmetric 2-sphere this is 

ds' = r'(d^' + sin2^d(/)'), (101) 

and thus y/deth — sin 6 giving the famihar 47rr^ area. Now if we ask how this area 
varies as we move from one slice of the foliation to another we are led to consider 

-f^ = ^V„^^ = f'VaAH = CtAH: (102) 

dA dA 

where A is a parameter labelling the foliations, and t"' is chosen to be tangential to the 
horizon, but normal to the foliation and always future-pointing. To relate this expression 
to the expansions of the ingoing and outgoing null rays rewrite the above as 

4-7r 

ATrCtrl = CtVdeth = 27rri h'^ jCthn- (103) 

sint' 

Now identifying the intrinsic metric with the projection tensor qab and writing Qab = 
&a &b + <Pa <Pb whcrc 9a and 0a are two orthonormal directions tangent to the 2 surface, 
we get 

^ = 4nrl q-'VM. (104) 

Here can always be written as a linear combination of and 

f'^Br + Cn" (105) 

where B and C are the coefficients and depend in general on A. However on a given 
slice, for t"- to be future-pointing, there are only five cases to consider: 
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t" outside r 


t"- oc t 


between i"- and n"- 


t"- oc 


inside n"" 


B>0 


B>0 


B>0 


5 = 


B <0 


C <0 


C = 


C>0 


C > 


C>0 



Additionally, on the horizon, 



dX 



^nriiCOn. 



(106) 



Therefore, since 9n is negative at the horizon by definition, the area Ah will be increasing 
in the first case, constant in the second case, and decreasing in the other three cases. 



9. Discussion 



We have demonstrated a simple framework capable of handling black holes with rapidly 
evolving horizons. The only physical simplifying assumption was one of spherical 
symmetry, and the two mathematical simplifying assumptions were the adoption of 
a spherically symmetric slicing of spacetime, and the use of a specific nonsingular 
coordinate chart (Painleve-Gullstrand coordinates). The framework does not assume 
asymptotic fiatness or stationarity, and we argue that these concepts are not required 
to define a useful concept of horizon. While exact spherical symmetry may not be 
an appropriate approximation to describe astrophysical black holes, we feel that this 
framework is useful in clarifying the basic concepts. In particular once someone 
has understood the simple nature of the evolving horizon in Painleve-Gullstrand 
coordinates, it is difficult to then become confused in other coordinate systems — the 
"frozen" nature of the horizon in Schwarzschild coordinates is then very clearly seen to 
be a simple coordinate artefact due to that coordinate system becoming degenerate on 
the horizon. 

Furthermore, this framework is useful in that any investigation of the black hole 
evaporation process must ultimately consider the essentially dynamic nature of the 
problem. While a particular choice of coordinates was used to facilitate some of the 
calculations, the simplicity of the results should in no way be seen as an exclusive 
property of the choice of coordinates. We have made no attempt to address the issue of 
singularity formation although Penrose's proof will remain valid within its assumptions. 
In this respect it is important to note that the current work assumes a classical spacetime 
everywhere, and one would expect it to arise as some sort of limiting approximation 
to the models of Ashtekar and Bojowald who consider the explicit breakdown of the 
spacetime manifold description. To avoid a singularity (as we certainly have trapped 
surfaces) we must either violate the NEC (or have closed timelike curves, or not have 
a manifold). Conservative approaches to avoiding singularities focus on the energy 
conditions and their violations. More radical approaches (either within loop quantum 
gravity or string models) effectively dispense with the spacetime manifold at sufficiently 
short distances. Still, before adopting the more radical approaches it is useful to see 
how far the standard manifold picture can be pushed. 
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One clear message that should be extracted form the current article is that event 
horizons [at least in their precise mathematical definition as absolute horizons] are not 
essential for doing interesting black hole physics. Event horizons are mathematical 
abstractions that in many ways encode and presuppose too many technical assumptions 
— assumptions which may or may not have anything to do with reality. For physical 
black holes, evolving either due to accretion or Hawking evaporation, other notions 
of horizon seem to be more appropriate. Apparent horizons are sometimes useful but 
are also known to possess several technical deficiencies. Isolated horizons are in many 
ways "too close" to being true event horizons (since no matter is allowed to cross 
the isolated horizon it cannot either grow or shrink). Hayward's trapping horizons, 
Ashtekar's dynamical horizons, and the evolving horizons of this article, all capture key 
aspects of back hole physics without presupposing the existence of an event horizon. We 
would argue that the key feature of the evolving horizons of this article is their relative 
simplicity and clarity, enabling one to build up a simple physical picture and a clear 
intuition. 



Appendix A. Coordinate basis 

For completeness, here are the coordinate basis components of the Einstein tensor: 

G« = ^ (^1 - — J + ^^f^r■ (A.l) 

2m 2m'c Phn 

Gtr^ — -—^\l—] (A.2) 



2m' 2c' 2m I , 

rc r'^c ^J2m/r 

Gee = G^^ r^; (A.4) 

G^4, = G^^ sin' 9. (A.5) 



At the evolving horizon 



Gtt^^] (A.6) 



2{m-m'c) _ 



rc r^c 



and 



c 

This then easily yields (on the evolving horizon) 



2m 



r^c 



which is a simple consistency check on the entire formalism. 
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Appendix B. Null decomposition 

Let the indices i,j, . . . only take on the values 0, 1 (so they lie in the t-r plane). Then 
spherical symmetry is enough to imply that 

Tij = To gij + T+ ii ij + T_ Uj. fB.l 

Then in terms of our natural tetrad 



-p + Pr ^ _ P + Pr - 2/ ^ _ P + Pr + 2/ 



In terms of the Einstein tensor 

2m' m 








G — + 

2r2c,y2^ ' 2rc 



At the apparent horizon 

2m' m 







9 ' 

m c' 



^ 2r^c rc 



B.2 
B.3 
B.4 
B.5 

B.6 
B.7 
B.8 



If wc now go to a general radially infalling tetrad = v + /3c, then Gq is independent 
of /? whereas the null vectors £ and n are Doppler shifted 



so that 



^ i±|; G_ ^ G_ i^. (B.IO 



In particular the product G_ is invariant under choice of radially infalling tetrad. 
This again demonstrates the while many statements that one might wish to make 
concerning the evolving horizon arc tetrad dependent, there are also a number of useful 
quantities that are tetrad-independent. In particular, note that if the evolving horizon 
is momentarily static, then since m = we have G- — 0, and the stress-energy on the 
horizon exhibits a type of "enhanced symmetry" [31] in that it takes on the form 

TiJ = Tog^,+T+eaJ■ (B.ll) 

Under the significantly stronger assumption that the evolving horizon is static, at least 
for some finite time interval, and has been static long enough for a bifurcation 2-surface 
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to form, then since £ — > on the bifurcation 2-surface, we have the much stronger result 
that 

Tij ^ To Qij (B.12) 

on the bifurcation 2-surface. This makes it clear that the "enhanced symmetries" of 
reference [31] are very much related to the existence of both a finite-time translation 
symmetry and a bifurcation 2-surface. The enhanced symmetry (B.12) is already 
reduced as one moves away from the bifurcation 2-surface, and even the reduced 
symmetry (B.ll) is at best only approximate once the horizon starts evolving. 



Appendix C. Special case: c{r,t) = 1 

There is a special case of the formalism in which technical computations simplify even 
further — that is if we set c(r, t) — 1. This is not a situation motivated by any particular 
physical principle: Although it is true that c(r, t) — 1 for both the Schwarzschild and 
Reissner-Nordstrom solutions, there does not appear to be any deep reason for this. 
(The analogous property also holds for the Doran form of the Kerr solution [32] , though 
again there seems to be no deep physical reason underlying this.) Nevertheless, the 
tremendous simplifications attendant on setting c{r,t) = 1 arc so significant, and do 
not seem to undermine the main points we are making regarding the non-perturbative 
evolution of horizons, that a brief summary may be warranted. (See the analogous 
discussion by Hayward in [14].) 

Suppose we consider the metric: 

ds^ = -[1 - v{r, t)^]dP + 2v{r, t) dr dt + dr^ + dQ^ (C.l) 



and set v{r) = \/2m{r, t)/r. Then the discussion of the evolving horizon is largely 
unaltered though now we have the simplification 

= = -v\r, t) = . (C.2) 

In fact, we now have 

Gahrt = GahU^n', (C.3) 
which is enough to tell us without detailed computation that 

GabV^S' = 0. (C.4) 

Indeed the "comoving" tetrad defined by V"^ and S"' is now geodesic, and the stress 
energy tensor in this tetrad basis simplifies to: 

^"8vr-47rr2' ^^'^^ 

/ = ^ = 0; (C.6) 
^ Grr ^ w! m 1 
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The only "complicated" component of the Einstein tensor is the transverse one: 



Gn 



Pt 



Stt 



Stt 



m 



Snr 



+ 



2m 



m 



?>2'Kvn?r 



\m — rm + 



1 



16m7r 



-m 



At the apparent horizon the only real simplifications are 



Pr,H = - 



m 



m 



H 



and 



Pt,H 



—mm" + mm' 



m 



;i-2m')/4' 



47rr^ 



(C.8) 



(C.9) 



(C.IO) 



H 



This particular class of spacetime geometries is slightly easier to deal with than those 
considered in the bulk of the article, and so is perhaps of some mathematical and 
pedagogical interest. However it should be noted that the more fundamental issues of 
interest in this article are the way in which we have explicitly shown how a suitable 
choice of coordinates completely side-steps the "frozen" nature of coordinate time at 
the Schwarzschild horizon, and allows one to non-perturbatively formulate questions 
about the production and decay of evolving horizons. 



Appendix D. Vaidya Spacetime 

The standard expression of the accreting Vaidya spacetime [5, 6] in advanced null 
coordinates is 

ds^ = - (l - + 2 d^; dr + dn\ (D.l) 

This spacetime represents collapsing null radiation, and the function m{v) is an arbitrary 
non-decresing function of its argument, f Note that the mass function is a function 
only of the advanced null coordinate v. % We can transform to Painleve-GuUstrand 
coordinates {grr — 1) by writing the null coordinate as a function of Painleve- 
GuUstrand coordinates t and r: that is, write v — v{t,r), whence 

d.= ^dt + ^dr. (D.2) 

at or 

f There is a related time-reversed "shining star" solution for outgoing null radiation, which can be 
written as 

ds^ = ~(l- ^HM] du^ -2dudr + r^ dQ^ 



where u is a retarded null coordinate and m{u) is a non-increasing function of its argument. This can 
also be used as a model for a rather specific class of evaoprating black hole. The comments below can 
immediately immediately be carried over to this case by suitable changes in notation. 
I Unfortunately the null coordinate v occurring here has nothing to do with the function v{r,t) 
previously used to discuss the Painlcvc GuUstrand form of the metric. The usage and notation is 
unfortunately standard, and care must be taken to keep the two concepts distinct. 
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2m{r{f.r))\ f dv 
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+ 



+ 



-2 1 



- 1- 



2m{v{i,r))\ / dv 
r J \dt 

2m{v{i,r))\ f dv\' 
\dr J 



dv\ f dv 
^ 1 +2( — 
or 



dv 
dr 



dt 



dt dr 



To obtain Painleve-GuUstrand coordinates we demand Qrr = 1 so 



- 1- 



2m{v{t, r))\ f dv\ 



\dr ) ~'~ ^ 



dv 
dr 



1. 



(D.3) 



(D.4) 



Solving this quadratic yields a first-order differential equation relating the coordinates 
{v, r) and (t, r): 



dv 
dr 



1± 



2'm(v{t ,r)) 



1 - 



2m{v{t,r)) 



It 



2m{v{i,r)) 



(D.5) 



Since v is by construction regular at the evolving horizon we must take ± — > — and 
^ — > + above. Putting this back in to equation (D.3) yields the Painleve-GuUstrand 
form of the Vaidya solution 



ds' 



- 1- 



2m{v{t, r)) 



dv 



dp 



+ 2 



2m{v{t, r)) f dv 
dt 



dt dr + dr^ + r^dfl^ 



(D.6) 



In the original null coordinates the evolving horizon occurs at 

rniv) = 2m{v), (D.7) 

and there is clearly only a single horizon for any value of the null coordinate v. In 
terms of the Painleve-GuUstrand coordinates {t, r) the location of the horizon is given 
implicitly by 

rH{v{i, r)) = 2m{v{i, r)), (D.8) 

where v{i, r) is constrained by the differential equation (D.5). If the horizon is spacelike 
this condition can in principle have many solutions for fixed t, as the horizon can in 
principle move back and forth in Painleve-GuUstrand time. Consider in particular the 
quantity dm{v(t, r))/ dr: 

dm dm dv dm 1 
dr dv dr dv , /2m(t,r) 



(D.9) 



Then as 2m — > r this becomes 

dm I 1 dm 
aT^""^" 2d^- 



(D.IO) 
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In order to violate equation (34) we would need m'^j > 1/2. But since in the Vaidya 
solution dm/dv can be any non- negative function, this condition can certainly be 
violated. In short, for a spacelike evolving horizon "outermost" in the sense of the 
{v,r) null coordinates may not always coincide with "outermost" in the sense of the 
{t, r) Painleve-GuUstrand coordinates. We emphasise that this is not an inconsistency 
in the formalism, merely one of the interesting coordinate artefacts one has to keep in 
mind. 
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